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Abstract 



We present a systematic expansion in the ratio between the level spacing 
and temperature and employ it to evaluate differences between statistical 
mechanics and thermodynamics in finite disordered systems. These differences 
are related to spectral correlations in those systems. They are fairly robust 
and are suppressed at temperatures much higher than the level spacing. 
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Typeset using REVT^ 



One of the major results that emerged from the work in the field of mesoscopic physics was 
the understanding that physical quantities may vary qualitatively depending on the averaging 
procedure employed, with a particular emphasis on differences between the canonical (CE) 
and the grandcanonical (GCE) ensembles. This observation became particularly apparent 
in the study of persistent currents in normal rings and has been extended to other 
thermodynamic and transport properties. A peculiarity of mesoscopic systems 

is the existence of an energy Ec associated with spectral correlations on scales larger than 
the mean level spacing, A. Earlier works [PHT^ have addressed differences between the 
ensembles at energies (or temperatures) comparable with A. Here we shall study differences 
between the CE and GCE which persist up to energy scales much larger than A. The term 
"canonical" deserves elaboration. One considers an ensemble of finite systems (typically 
conducting), which are macro scopically similar but differ from each other by their impurity 
configurations, details of the boundaries etc. We now assign A^*^*) electrons to the i*'* member 
of the ensemble; A^^*^ is selected arbitrarily, and is uncorrelated with the energy spectrum of 
the system. Moreover, it is unchanged as we vary some external parameter, x (representing 
e.g., an Aharonov-Bohm flux). In particular one may select A^*^*-* = N = const for all i. 

The idea introduced in Refs. 0] was to represent each member of the canonical ensemble 
as a grandcanonical system with an effective chemical potential, fx^^^x), which is sample 
specific and depends on the parameter x. A canonically averaged quantity can be obtained by 
expanding around a "wrong" , average, chemical potential, /i (which is independent of i and 
x), and then calculating the expansion terms (taken at /i, i.e., performed grandcanonically) . 
This procedure, leading to remarkable differences between CE and GCE averages, is clearly 
warranted at zero temperature, T=0. In that case differences between the CE and the GCE 
are solely due to averaging over quenched disorder: each system individually can be equally 
described by either A^ or fi^^\x), such that ||13| 



N{fx\x))=N. (1) 



Here n^'\n, x) is the (sample specific) thermodynamic (TD) potential. The procedure al- 



luded to above has also been employed at T > 0, defining the sample specific through 
Eq. (|1]). In that case a canonical system characterized by a fixed number of particles, A^, is 
replaced by a grandcanonical system whose expectation value of the number of electrons is 
(at any given value of x). According to this picture differences between the CE and the 
GCE are again due to quenched disorder (but may depend on temperature). The starting 
point here is the TD relation 

F(^) ( A^, x) = fi(^) (/i, x) + fiN (2) 



(F*^*) is the free energy), leading to the identity: 



N dx 



dx 

where A^ is related to fx^^^x) through Eq. (1). We should recall, though, that the TD rela- 
tions, Eqs. d^) and (§), are only approximate (see Eq. (||) below) when it comes to statistical 
mechanics (SM). It is the latter which should be employed to calculate expectation values 
of equilibrium observables. F*^*^ should be derived microscopically, and not by calculating 
^(0 microscopically and then employing Eqs. (2), (^. To understand what is missing in 
the procedure described above, one may take as an example a system of non-interacting 
electrons: the occupation probability for a canonical system is not given by the Fermi-Dirac 
function with an effective fi Differences between canonically and grandcanon- 

ically averaged quantities should reflect two elements: (i) sample-to-sample fluctuations 
due to quenched disorder; (ii) the deviation of SM from TD resulting from the fact that 
the probability to find the system in a given quantum state differs between the GCE and 
the CE. The latter element, manifested in corrections to Eqs. (^ and (0), occurs in finite 



systems and only at finite T. 

The purpose of the present work is to elucidate some basic questions concerning the 
thermodynamics of finite systems. We report on a systematic study of differences between 
the CE and the GCE, yielding both the contributions due to quenched disorder and due to 
differences between SM and TD. We specifically consider non-interacting electrons moving 



in a diffusive disorder at T > A. We demonstrate our approach by evaluating both the 
persistent current and the heat capacity. We indicate how both types of contributions reflect 
spectral correlations in the system. Our analysis also points out how the thermodynamic 
limit (where differences between SM and TD vanish) is approached. Our general scheme 
should allow for the study of other regimes of disorder (e.g. the dirty ballistic), other types 
of averaging (e.g. over energy), and the inclusion of electron-electron interactions. 
The canonical partition function is given by 

Z»(iV)=Tr{5(iV-iV)e-^^*''} = ^ J d/i e-^(^'''(^)+^^) , (4) 



-ittT 



where hatted quantities are operators, e~^^ = Tr{e^'^*-^^^^^} and P = l/fc^T. Since the 
eigenvalues of N are integers, the 5-function is in fact a Kronecker function. We note that 
evaluation of the integral in Eq. (^) within saddle point approximation |T5|JT6[] , leads to the 
TD relation Eq. (0). The saddle point JiS'^\x) is sample specific. Instead we expand around 
a constant /I, defined by {N^'^\fl)) = N, where (. . .) denotes averaging over realizations of 
the disorder. Defining fT^ /i - /i = iVATr and V^^ = {n\T)-^{-M'Yl'^d''-^5N'^^ / dpL^'^^ 



we obtain 



2ttT J 



.n=l 



(5) 



Here 5N^^\jj,, x) = —dVl^'^ydfi — N is the sample specific deviation from the mean number 
of electrons. For A/T ^ 1 we may replace the integration limits in Eq. (|^) by ±oo. 
Expanding the term exp [X]„ V^^* V"] , we obtain a diagrammatic expansion where the {V^^*-*} 
play the role of n — th order vertices (see below). We stress that the expansion around ft is 
not a standard one: ft does not represent the sample specific saddle point. This expansion 
is justified a-posteriori (see below). Employing linked cluster diagrammatic expansion []TS|, 
the free energy may be written 

F^'\N) = {n^'\fi) +fiN) + - In f ^ j - {all connected diagramsj^^^ . (6) 



The first term in Eq. yields sample specific grandcanonical observables. The second term 
due to Gaussian fluctuations around the expansion point at fl, is a T-dependent sample- 
independent contribution. 

We next consider the diagrammatic expansion of F»(Ar). This expansion involves inte- 
grals over T. In a somewhat artificial, but convenient, analogy with standard perturbation 
theory, we refer to the term (27r)^^/^ /rfr exp{— r^/2}r^ = 1 as a contraction or a "free 
propagator", symbolically (tt). Similarly a factor (27r)^^/^ /(irexp{— r^/2}r^ = 3 involves 
3 different ways of contraction, (tttt) = 3. Our "propagators", or statistical lines, carry 
neither energy nor momentum, and are employed to identify contributing diagrams and to 
evaluate combinational factors. Each propagator is then represented by a "statistical line" 
and each factor Vn corresponds to an n — th order vertex, representing an electron loop 
with n scalar vertices. Examples of connected diagrams contributing to F^^\N) are shown 
in Figs. ID, ^ These are skeleton diagrams which are yet to be dressed by impurity lines, 
interaction lines, etc. 

The following set of rules applies for the calculation of the skeleton diagrams: (1) consider 

a connected diagram consisting of p loops with rii, . . . ,np scalar vertices respectively. We 

require = 2k, with k being the number of statistical lines, k > p — 1. (2) Each 

statistical line carries a factor (—A). (3) Each n- vertex loop carries a factor ^ q-^-i ■ 

P- 

(4) Each diagram carries a factor T^, where w = k~p+l>0. (5) A factor of 1/m! should 
be assigned to any subset of m loops having the same number of vertices. (6) Each diagram 
should be multiplied by a combinational factor reflecting the number of different ways to 
interconnect vertices by lines, which result in the very same diagram. 
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FIG. 1. The w = family of skeleton canonical diagrams with k = 1,2,3. Zigzag lines are 
"statistical lines"; full lines - electron propagators; black dots - scalar vertices. The first diagram 
upon averaging over disorder yields Altshuler-Shklovskii result. 



We now divide the skeleton diagrams into families, characterized by the index w. Exam- 
ples are depicted in Figs. 0, Ij. It is understood that disorder averaging has to be carried out 
subsequently. Let us first consider the w = family. Fig. 0. It turns out that these are the 
contributions (describing differences between the CE and the GCE), obtained when Eqs. (0) 
and (0) are assumed to hold and each member of the canonical ensemble is assigned an 
effective sample specific chemical potential: the w = family represents contributions due 
to quenched disorder, but not due to differences between SM and TD. (It is the only con- 
tribution which survives at T = 0). After some algebra one obtains (hereafter we suppress 
index (i) and consider only ensemble averaged quantities) 

Exactly the same result may be obtained by the direct solution of Eqs. (|l]) and @) [|T9| . 
The first {k = 1) term in the sum corresponds to the two loop diagram in Fig. |l|. Upon 



averaging over (diffusive) disorder, it yields the Altshuler-Shklovskii term pO|, which has 
been employed in Ref. 0. The k > 2 terms are given by complete derivatives with respect 
to yU, and are negligible upon averaging (being small in the parameter A//i). This provides 
us with an a-posteriori justification of the diagrammatic expansion: sample specific terms 
in this expansion are not necessarily small, but the ensemble average is well-behaved. 

We next include the w > 1 families too. These yield contributions due to deviations of 
SM from TD. It can be shown that for the regime where we employ our expansion, T > A, 
the leading term (in A/T) of each family is represented by a two-loop diagrams (cf. Figs. 
|l], 1). Evaluat ion of these diagrams (neglecting full derivatives with respect to fi) results in 

_ A - (AT)- /(d-6Ny\ 



The w = term is the Altshuler-Shklovskii quenched disorder contribution ||2^; w=l is the 



leading SM vs. TD contribution 



[21 
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FIG. 2. Skeleton diagrams of the w = 1 family (up to k = 3). 
To evaluate these terms we define the correlator K and its Fourier transform K: 



^ +00 



K{e - e') ^ A' - ^- J dtk{t)e^'^^-^'^l^ , (9) 

— oo 

where v{e) is the sample specific density of states and t is the dimensionless time (in units 
oih/ A). We can write 



d'^SNW d"" d"" rr dede' 



dp,"" J I dfi"" dfi"" J J A2 



(10) 



fi'=ii 



where / is the Fermi-Dirac function. Changing variables to ^ = e — e', 77 = (e + e')/2, 
performing the integral over 77, and Fourier transforming with respect to ^, we obtain 



(«..-,o^) = ^i:t-^ M M-tcta- *^ • (11) 



Eqs. (^ and (11) form the basis for the analysis of the various corrections to the GCE 
averages, and depict the dependence of these corrections on spectral correlations. We study 
two examples comparing the w=0 (quenched disorder contribution) and w=l (leading con- 
tribution due to differences between SM and TD). 

The average persistent current in the canonical ensemble, {I)ce, is obtained by deriving 
(F) with respect to the Aharonov-Bohm fiux, 0. The fiux dependent part of the time 
correlator for quasi-one-dimensional rings is given by 



t / 1 ^ 
^('^) = \ -A — m E cos47rp0/0o , (12) 



IT y 47r(7|t| 



p=i 



where g = Ec/A ^ 1 is the dimensionless conductance {Ec is the Thouless correlation 
energy) and (pQ = hc/e. Only even harmonics appear in Eq. (|12|). Expanding {I)ce = 
-^psin47rp0/0o, we find that the quenched disorder contribution is 



A 



2 

A < T < Ec/p^ 



% I /^e"V ^ T > Ejp' 



E. 



c 



while the leading SM vs. TD correction is given by 



jSM-TD _ ^1 

^ <Pog 



15v^C(5/2) , ^ ^ , , 

1 T - /^r) 



The dependence of the SM-TD contribution on temperature and fiux is significantly different 



from the contribution due to quenched disorder (Eq. (|13D ), which has been found previously. 
The former shows that differences between SM and TD decrease slowly (algebraically) with 
temperature up to T ~ ^ A, and then are suppressed exponentially. 

Average canonical heat capacity {{C)ce)- It is given by {C)ce = — /3^<9^(/3(-F))/c?/3^ and 
is written as a sum of contributions (cf. Eq. (|^)): 

{C)cE = {C)gce + + {SCn + (5C^*^-™) . (15) 

The first term is the grandcanonical contribution, which for a degenerate gas of non- 
interacting electrons is The term due to Gaussian fiuctuations yields — |. This contri- 
bution can be reinterpreted as a shift of the grandcanonical temperature (T > A) towards 
a lower temperature, T T — forced by the canonical constraints. A similar shift has 
been found in Ref. 0. The next terms in Eq. (pIS]) are evaluated employing Eq. (0). We are 
interested in energies larger than A, hence t < 1. The time correlator has two interesting 
regimes (random matrix theory and Altshuler-Shklovskii respectively): 

I/I f 1 g-^ <t<l 

Kit) = ^ , (16) 

i {A7ig\t\)-'^/^ ATei <t<g-^ 

where h = 1,2,4 for the orthogonal, unitary and symplectic ensembles respectively. Here 
Tel is the elastic mean free time. Eq. (|r^) does not account for the crossover regimes. This 
leads to 
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1 A I ^ A < T < 



and in the leading order in A/T ||22| 

f 3C(3) 



- — A<T < Ec 

in 



rp.d/2 (18) 



r/. ( - ) T > 



c 



In this example the SM-TD contribution is parametrically smaller than the contribution due 
to quenched disorder. Here too, differences between SM and TD decay slowly over scales 
larger than A: as a power law, at least over energy up to h/Tei- The results depicted in 
Eqs. (p!?!) and (|l^), as well as an evaluation of {6C^^^~'^^) for a Poissonian spectrum 
lead us to conclude that the less rigid the spectrum the slower the decay with temperature. 

In summary, we have presented a systematic expansion in A/T describing differences 
between statistical mechanics and thermodynamics of finite systems. These differences, 
related to spectral correlations, are found to be fairly robust and slowly suppressed with 
T. The thermodynamic limit where they are suppressed is attained at temperatures much 
higher than the level spacing. 
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Science Foundation (BSF) and the Israel Academy of Sciences. 
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